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Abstract 

We examine resonant slepton production at the LHC with gravitinos in the final 
state. The slepton undergoes gauge decay into a neutralino and a lepton, the neu- 
tralino decays into a photon and a gravitino. By measuring the transverse masses 
of the 7G and the Z7G subsystems it is possible to accurately reconstruct both the 
slepton and the neutralino masses. In some regions of parameter space the slep- 
ton decays directly into a lepton and gravitino, giving an identical experimental 
topology to W production (l$ T ). We present the novel matrix element squared for 
lepton-gravitino production. A peak in the tail of the lepton-missing momentum 
transverse mass distribution of the W provides a signature for the process and an 
accurate measurement of the slepton mass. We display the search reach for the LHC 
and 300 fb _1 of integrated luminosity. 



1 Introduction 



This letter is devoted to the study of the signals at the Large Hadron Collider 
(LHC) due to a supersymmetric generalisation of the Standard Model (SM) 
which (a) violates i?-parity, and (b) has an ultra-light gravitino in its spectrum. 
Despite the negative results that experimental searches for supersymmetric 
particles have yielded so far, the philosophy underlying these searches has been 
rather exclusive. The so-called minimal supergravity inspired .R-conserving 
model of minimal supersymmetry or the constrained minimal supersymmetric 
standard model (CMSSM), as it is sometimes called, has been elevated to 
the status of a paradigm in the quest for supersymmetry at colliders. The 
adherence to this model to the exclusion of a myriad of other possibilities is 
not desirable. This is motivation enough to consider other models which relax 
some of the assumptions made in the CMSSM. 
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The quest to depart from the narrow confines of the CMSSM inexorably leads 
one to question the assumption of -R-parity conservation as well as seeking 
alternatives to the mass spectrum of supersymmetric (SUSY) particles pre- 
dicted by the minimal supergravity models. Of these two paths, the second 
is certainly more interesting, connected as it is with the fundamental issue 
of SUSY breaking. Of particular interest is the question of the lightest su- 
perparticle (LSP) because the decay patterns of the heavier particles into the 
LSP decide several of the collider search strategies for supersymmetry. In the 
preferred models, the LSP (most often the neutralino) is assumed to be in the 
mass range just above the reach of present experiments i.e. of 0(100) GeV 
in mass. All the other sparticles, are by definition, heavier than the LSP (the 
next-to lightest we denote the NLSP). A dramatic alternative is to have an 
ultralight gravitino as the LSP. The gravitino, in some models of supersym- 
metry breaking such as gauge-mediated supersymmetry breaking [1], can be 
as light as 10~ 3 eV: other sparticles then decay into the gravitino and this 
can significantly alter the strategies for supersymmetric particle searches at 
colliders. The other scenario wherein the LSP is not stable is that of -R-parity 
violation where the decays of the neutralino yield final-states with jets, leptons 
or missing- E T . But -R-parity violation has another interesting implication, ie 
the sparticles can be produced singly. If this single sparticle production pro- 
cess proceeds via the s-channel, then the cross-section is resonantly enhanced 
and may be significant even for small values of the -R-parity violating coupling. 

An anomaly in the CDF experiment in the production rate of lepton-photon- 
$ T in pp collisions at y/s = 1.8 TeV was observed using 86.34 pb _1 of Tevatron 
1994-95 data [2]. While the number of events expected from the SM is 7.6 ± 
0.7, the experimentally measured number corresponded to 16. Moreover, 11 
of these events involved muons (with 4.2 ± 0.5 expected) and 5 electrons 
(with 3.4 ± 0.3 expected). These anomalous events can be simply understood 
[3] in terms of a supersymmetric model with the following features: (1) the 
model is -R-parity violating with an L-violating A' 211 coupling, and (2) the 
supersymmetric spectrum includes an ultra-light gravitino of mass ~ 10~ 3 
eV. The resonant production of a smuon via the .R-violating coupling, its 
decay into neutralino and a muon and, finally, the decay of the neutralino 
into a gravitino and a photon leads to the /ry-^ T final state studied in the 
CDF experiment. The range of smuon and neutralino masses relevant to the 
explanation of these anomalous observations of the CDF experiment is such 
that most of this range will be explored at the Run II of the Tevatron. In the 
event that this signal is not seen at Run II it will rule out the model at the 
lower end of the neutralino and smuon masses. 

For heavier smuon and neutralino masses (above 250 GeV, roughly), the afore- 
mentioned Run I signal would be a statistical fluke and will probably disappear 
in Run II data. In that case, the LHC can be expected to discover and measure 
the sparticle masses. In this paper, we perform a study of the ability of the 
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LHC to perform these two tasks, identifying the sensitive observables. 



2 The model 

We assume throughout this paper that of all the i?-violating couplings, only 
one lepton-number violating coupling is significantly large while the others are 
too small to have observable consequences. The coupling appears as 

Kjk L iQj D k C W (1) 

in the superpotential (where k are family indices, gauge indices have been 
suppressed and Lj, Qj, denote the left-handed lepton doublet, the left- 
handed quark doublet and the charge-conjugated right-handed down quark 
superfields respectively). This coupling allows for the production of a slepton 
from the initial state of quarks. The slepton thus produced could decay again 
via the i?-violating coupling into two quarks. However, if the R- violating cou- 
pling is small, the existence of an ultralight gravitino of a mass ~ O(10~ 3 ) eV 
drastically alters the decay mode of the slepton. The slepton overwhelmingly 
decays into a lepton and a (bino-dominated) neutralino, with the latter decay- 
ing into a photon and a gravitino resulting in a lj$ T final-state. The Feynman 
diagram for the process is shown in Fig. 1. Since < mj L we should also 
expect signals from sneutrino production. In the case of sneutrino production, 
the process qq' — > v — > 7GV gives rise to a 7^ T final-state. We will not study 
this final state in this paper, since the background depends crucially on cos- 
mic ray events which are difficult to estimate. In our model, the only other 
light sparticles are the neutralino (which is lighter than the slepton) and the 
ultra-light gravitino (which is as light as 10 -3 eV). We enforce degeneracy 
between the first two generations in order to avoid flavour changing neutral 
currents. Other sparticles do not play a role in this analysis, and are set to be 
arbitrarily heavy. 

Such a light gravitino materialises naturally in models of gauge-mediated su- 
persymmetry breaking (GMSB)[1]. However, in minimal GMSB models the 
chargino and the second-lightest neutralino are not much heavier than the neu- 
tralino. This is certainly not desirable for our considerations because it predicts 
large jets+7+^ T rates which are not seen by experiments 1 . The question that 
naturally follows is how we can deviate from the minimal versions of GMSB 
models and obtain the pattern of supersymmetric masses that we require. It is 
interesting to note that such a mass spectrum can arise in an alternate model 
of GMSB which is obtained from compactifying 11-dimensional M-theory on a 

1 This observation has been made earlier in the literature [4] in the context of the 
GMSB-based explanation [5] of the eejj$ T CDF event [6]. 
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Fig. 1. Feynman diagram of resonant slepton production followed by neutralino 
decay. 

7- manifold of G2 holonomy [7]. However, the detailed low-energy predictions 
of this model have not been worked out so for our analysis we simply work 
with the low-energy model described above. 



3 Simulation Results 

For our study of the process shown in Fig. 1 at the LHC (pp collisions at the 
y^s = 14 TeV), we have chosen to work with the following default set of model 
parameters (unless indicated otherwise): 

• Gravitino mass, = 10~ 3 eV, 

• R- violating coupling A' = A^u = 0.01, 

• tan/3 = 10, 

• sparticle masses (m x? ,m r ) = (120 GeV, 200 GeV) or (200 GeV, 500 GeV) 
GeV ("low mass" and "high mass" scenarios) respectively. 

The choice of using X' 2n rather than some other flavour combination is ar- 
bitrary. Our results can be easily generalised to other i?-violating couplings 
and we will comment upon this below. Note that the chosen value of A' 211 
is much smaller than the bound coming from R n = T(tt — > eu)/(ir — > fiu) 
[8]: A 211 < 0.059 x ^ [9] even for a squark mass of 100 GeV. However, 
since the squark mass is arbitrarily large in our model, this constraint is not 
relevant to our analysis. The .R-violating decay of the slepton is possible but 
constrained, in principle, by the Tevatron di-jet data [10] which exclude a 
a.B > 1.3 x 10 4 pb at 95% C.L. for a resonance mass of 200 GeV. However, 
in practice this does not provide a restrictive bound upon our scenario as long 
as the i?-violating coupling is sufficiently small, < 0(1). Moreover, the di-jet 
bound is not restrictive because it suffers from a huge QCD background. By 
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restricting A^n and the gravitino mass to be small, we avoid significant rates 
for the possible R- violating decays of Xi ~^ M3 or Xi ~^ u 33- Xi ~^ lG is the 
dominant channel. 

We use the ISASUSY part of the ISAJET7.58 package [11] to generate the spec- 
trum, branching ratios and decays of the sparticles. For an example of param- 
eters, we choose (in the notation used by ref. [11]) tan/5 = 10 and A t)T> b = 0. 
\i together with the flavour-diagonal soft supersymmetry breaking parameters 
are set to be very large. We emphasise that this is a representative point in 
the super symmetric parameter space and not a special choice. 

3.1 Slepton decays to lepton, photon and gravitino 

The signal has been simulated using HERWIG6.4 [12]. The SM background 
to the [i^E T and ^E T channels have been simulated using PYTHIA [13]. Be- 
cause the background took a restrictively long time to calculate, 10 times less 
background luminosity was simulated than signal luminosity. The background 
has then been scaled up by a factor of 10 for all results presented here, but it 
should be borne in mind that the statistical fluctuations within it are bigger 
than will be expected. For our simulations of the signal and the background we 
have used only selection cuts of E T > 25 GeV on the transverse energies of the 
muon and the photon. The same cut of 25 GeV is also used for E T . We have 
used the following cuts on the rapidity of the photon and the muon: |t7 7j/1 | < 3. 
There is an isolation cut between the photon and other hard objects o in the 



event of y (r/ 7 — r] ) 2 + (</> 7 — <p ) 2 > 0.7. Since the signal is hadronically quiet, 
we veto events with jets reconstructed with E T > 30 GeV and rjj < 4. Initial 
and final state radiation effects, as well as fragmentation effects are included 
in the background simulation. 

The E T distributions of 7, /i and E T in the simulated events are shown in 
Figs. 2 (a)-(c). It would be very difficult to claim a signal for the low or high 
mass scenarios based on these distributions, given uncertainties involved in 
the background calculation. However, the figures illustrate the fact that the 
signal to background ratio could be improved by tightening the Et cuts, to 50 
GeV for example. We leave the optimisation of cuts to a future, more detailed 
study. 

The transverse mass distributions of several particles i 



should possess a sharp peak at the mass of the particle which decays to i daugh- 





(2) 
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Fig. 2. Et distributions of (a) muons, (b) photons and (c) missing-E^ for h^/^t 
events that pass the cuts. 300 fb _1 integrated luminosity at the LHC is assumed. The 
purple (lighter) histograms display W^f SM background, the red (darker) histograms 
show the signal plus background for (m x o,mj) = (120 GeV, 200 GeV), whereas 
the blue (dotted) histograms display the signal plus background distributions for 
(m x o,mj) = (200 GeV, 500 GeV). 



ters. The Mt{[j>j$ t ) and Mt(7-£/ t ) distributions are displayed in Figs. 3a,b for 
the simulated high and low mass points and simulated SM background. 
In Fig. 3a, sharp peaks in the M T (/z7^ T ) are clearly visible at values of the 
smuon mass and will be detected above the SM W'-f background. The smuon 
mass could be accurately measured in this manner. The background gives a 
negligible number of events above M^(/i7^ T ) > 450 GeV. Fig. 3b shows that 
the signal peaks in Mt(j$t) (predicted to be at the neutralino mass) should 
be able to provide a measurement of the neutralino mass. 

In order to calculate the search reach, we use the signal S in the 4 highest 
peak bins (covering 20 GeV) of M T (/i7^ T ). The background distribution in 
these four bins is estimated by fitting a simple function to M T {^n$ T ) between 
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Fig. 3. Mr distributions of (a) /U7^ T , (b) for 300 fb _1 integrated 

luminosity at the LHC. The purple (lighter) histograms display W7 SM 
background, the red (darker) histograms show signal plus background for 
(m x o,m[) = (120 GeV, 200 GeV), whereas the blue (dotted) histograms display 
the signal plus background distributions for (m x o,mj) = (200 GeV, 500 GeV). In 
(a) , the dashed black line displays a log- linear fit to the background distribution for 
Mt = 150 — 400 GeV. In (b), the insert shows a linear scale magnification of an 
area of the plot. 




V 2U m- (GeV) 



Fig. 4. (a) Search reach for the ^7^t signal (as defined in the text) for 300 fb _1 
integrated luminosity at the LHC. (b) Relative production cross-sections for the 
different flavour choices of X' ijk couplings as a function of the slepton mass. 



150-400 GeV in Fig. 3a. We use 



B = 4exp [aM r (/i7$ T ) + b] 



(3) 



with purely \/l3 statistical errors, leading to a = —0.018 ± 0.001, b = 9.25 ± 
0.22. B is displayed in Fig. 3a as the dashed black line. We define the sensitive 
region of parameter space to be any region for which 2 S/ \f~B > 5 and S > 10. 
For 300 fb _1 , the search reach is shown as a function of smuon mass and Im- 
parity conserving coupling in Fig. 4a for the case of a neutralino NLSP. 

If a different R-parity violating (RPV) flavour coupling than X' 2U was used, 
the analysis would be much the same provided the coupling did not involve 
taus or top quarks. The main difference would be that, for a given size of RPV 
coupling, the production cross-section would decrease because one no longer 
probes valence partons in the proton. Fig. 4b shows that the suppression of 
the cross-section also depends upon the slepton mass, since different regions 
of x are being probed in the parton distributions. For a fixed slepton mass, 
the ratio of production cross sections cr(\ijk = c)/a(Xm = c) is shown. The 
parton distributions of top quarks are not known, so we do not include the 
possibility that the RPV coupling involves the top (A^). The figure shows 
that the production cross-section can decrease to less than 10% of its original 
value by changing the flavour structure of the RPV coupling to include heavy 
quarks. The number of produced events for different flavours of RPV coupling 
other than A' 211 can therefore be determined by multiplying the ratio of cross 
sections in Fig. 4b with the number of events in previous plots. One must 
then ask how the number of measured signal events depends upon the lepton 
flavour of the RPV coupling used. Efficiencies and backgrounds for the e and 
/i case will differ somewhat depending upon the experiment, and so changing 
the lepton flavour involved in the RPV coupling could quantitatively affect 
the result, r reconstruction efficiencies are significantly different to those of 
e, jj, and so there would be a large change in number of measured events in 
that case. 

There is no particular need to rely on a neutralino NLSP. Indeed, recent 
studies have found sleptons to be lighter than neutralinos in, for example, 
large regions of RPV mSUGRA space [14]. Later (for a different final state), 
we will consider the possibility of a smuon NLSP, and contrast it with the 
neutralino NLSP case. 



3.2 Slepton decays to gravitino and lepton 

We now turn to the decay I — > Gl. We ignore sneutrino production in this case 
because it would lead to an invisible final state. Neglecting lepton masses, the 
partial widths for the decays of I [1,8] are 

2 The statistical uncertainties on fitted a and b parameters make a negligible dif- 
ference to the final numerical results. 
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Fig. 5. Feynman diagram of resonant slepton production followed by decay into a 
lepton and gravitino 
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where b = eNu + gNi2(l/2 — sm 2 W )/ cos9 w and Ny denote neutralino mixing 
matrix elements [15]. T(l — > Gl) thus dominates for large mj and small rug. 



The Feynman diagram of the process in question is displayed in Fig. 5 and we 
calculate the spin and colour averaged matrix element squared to be: 



|M| 



2 _ \' 2 m\ 



36Mj,m}, (s - m r ) 2 + T 2 m 



2 ■ 



(5) 



where s is the centre-of-mass energy of the colliding qq' system, Mp = 2.4x 10 18 
GeV is the reduced Planck mass and T is the total decay width of the slepton. 

This new process is not simulated in HERWIG, so we perform the phase space 
integration by the method of Monte-Carlo integration ourselves. For this rea- 
son, the simulation was not of generated events and is therefore not subject 
to associated statistical errors. We have use cuts on the muon and missing 
transverse energy identical to the 7/^ T analysis, i.e. Ej. > 25 GeV and 

The M T (fiE T ) distribution is displayed in Fig. 6a for the SM W background 
plus signal in two cases: = 1.2 TeV, 1.5 TeV. If we use the default value 
of A' 211 = 0.01 and the neutralino is the NLSP, we obtain the smaller peaks 
which would not be prominent enough to claim a 5a discovery. If the smuon 
is the NLSP, the decay mode to neutralinos is not open, leading to a higher 
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Fig. 6. (a) Mt distribution of the \iG final state for two different values of the 
smuon mass and 300 fb _1 integrated luminosity at the LHC. The left(right)-most 
three peaks are for smuon masses of 1.2 and 1.5 TeV respectively, (b) Search reach 
(as defined in the text) for the \iG final state, a neutralino NLSP and 300 fb _1 
integrated luminosity at the LHC and A^n = 0.1. The search reach is contained to 
the left of the curve. 



branching ratio for jl — > G\x. The magenta (middle height) curves in Fig. 6a 
are evidence for this, and should be observable even for the small value of 
A211 = 0.01. The signal cross section scales as A^q 2 , and for A^n = 0.1 we 
obtain the one-hundred times larger peaks, either of which would be easy to 
detect on top of the SM W background, even for a neutralino NLSP. 

In this subsection, we define the search reach by the criteria that for S signal 
and B background events, S/ \fB > 5 and S > 10 in the 5 GeV signal peak 
bin of Mf(/i^ T ). These constraints lead to the expected search reach for 300 
fb _1 luminosity shown in Fig. 6b for a neutralino NLSP. For smuon masses 
that are too low, the huge background swamps the signal in the peak bin. For 
smuon masses that are too high, the signal production cross-section is too low. 
The search region is shown as the region to the left of the curve in Fig. 6b. 

As explained above, it is easier to detect the [iG final state if the smuon is the 
NLSP, because the relevant branching ratio is larger. Thus we expect a larger 
region of sensitivity. The search reach in two different planes of parameter 
space are displayed in Fig. 7 for the case of a smuon NLSP. For X' 2U = 0.1, 
the region of sensitivity is modestly larger than the neutralino NLSP case, as 
can be seen by comparing the blue curve in Fig. 7 with the one in Fig. 6b. 
However, there is no possibility searching for a peak with A' 211 = 0.01 with a 
neutralino NLSP, whereas a smuon NLSP is still detectable up to m 1500 
GeV, as displayed in Fig. 7a. 
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Fig. 7. Region of sensitivity to the Mt{^G) peak for 300 fb _1 at the LHC for a 
smuon NLSP. The region of sensitivity is (a) to the left hand side of the curves, (b) 
to the right hand side of the curves. 

4 Conclusions 



We have provided a basic first study of the search for resonant slepton produc- 
tion at the LHC with ultra light gravitinos in the final state 3 . We have not 
performed a detailed detector simulation, with associated mis-tagged back- 
grounds, Ifirp resolution and detection efficiencies. We started with a simple 
and loose set of cuts, and provided distributions which indicate how signal to 
background ratio can be improved. Specific cases of RPV coupling and other 
parameters are taken in order to be concrete, but the analysis applies over a 
wide range of parameter space. 

The current paper shows which are the sensitive variables (various Mt dis- 
tributions are the most useful) and displays the rough search reach one can 
expect. Luckily, signal peaks in My distributions mean that one does not 
have to know the backgrounds very well. The backgrounds can be measured 
and fitted away from the peaks, whereas the position of the peak provides 
an accurate measurement of the slepton mass. Although detector efficiencies 
and non-physics backgrounds will erode the search reach somewhat, it could 
be possible to flavour subtract some backgrounds. SM background processes 
from W — > Iv are lepton-flavour universal. On the other hand, if one RPV cou- 
pling is dominant over others, the signal only contributes to the production 
of one lepton flavour. Thus, subtracting the number of electron-tagged events 
from the number of muon-tagged events would greatly improve the signal to 
background ratio. The only pitfall of this approach is that it is possible that 

3 We note that resonant slepton production at the LHC has been studied before 
without decays into gravitinos in Ref. [16], for example. 
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the parameter space is such that one produces similar numbers of selectrons 
and smuons and therefore one could inadvertently cancel the signal. Although 
this only happens in a very limited region of parameter space, one cannot 
a priori rule this pessimistic scenario out and so we have not included it in 
our analysis. The process qq' — > I — > Gl has not been studied before, and we 
presented the matrix element squared together with the search reach and Mt 
distributions. 

We have concentrated on LHC studies, but of course the same ft — ► [iG pro- 
duction analysis could be applied to the Tevatron, with lighter smuons than 
those considered here. The important indicator would still be a peak in the 
Mt(h) tail, above M\y. Instead of smuon production, we could instead have 
considered squark production through a baryon-number violating coupling 
l/2\" jk UiDjD k in the superpotential. We would then have signals of j^$ T in 
the neutralino NLSP case 4 and j$ T ("monojet" signature). The sensitivity 
of the Tevatron and LHC to these processes is presumably less than in the 
slepton-production case, but remains to be calculated. In that case, the same 
Mt variables will be appropriate for detection, but some of the cuts may need 
to be harder in order to beat down a larger background. 
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